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ABSTRACT 



We analyze variations in the transit lightcurves of CoRoT-2b, a massive hot Jupiter orbiting a highly active G star. We use one transit 
lightcurve to eclipse-map a photospheric spot occulted by the planet. In this case study we determine the size and longitude of the 
eclipsed portion of the starspot and systematically study the corresponding uncertainties. We determine a spot radius between 4.5° 
and 10.5° on the stellar surface and the spot longitude with a precision of about ±1 degree. Given the well-known transit geometry of 
the CoRoT-2 system, this implies a reliable detection of spots on latitudes typically covered by sunspots; also regarding its size the 
modelled spot is comparable to large spot groups on the Sun. We discuss the future potential of eclipse mapping by planetary transits 
for the high-resolution analysis of stellar surface features. 
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1. Introduction 

The atmosphere of the Sun shows inhomogeneities down to the 
smallest scales currently accessib le to solar observati ons which 
are of the order of 50 km (e.g. IScharmer et alj |2002). Such in- 
tricate fine structure can also be expected in the atmospheres of 
other active stars. However, the best currently available stellar 
observations only resolve surface features down to the size of a 
few degrees on the surface, corresponding to several 10000 km 
on a main sequence star. 

The increasing number of known transiting extrasolar plan- 
ets offers an outstanding opportunity to study surface inhomo- 
geneities of their host stars with an unprecedented surface res- 
olution. As an example of this new technique, we present the 
analysis of one transit lightcurve of the plan etary system CoRoT - 

2, recently detected by the CoRoT satellite dRouan et al.Ul998l) . 
Our study indicates that under favourable conditions the CoRoT 
lightcurves allow the study of e.g. starspots and/or faculae down 
to a sub-degree scale on the stellar surface. 

Deformations of planetary transit lightcurves, attributable 
to dark spots , have been observed for several systems: 
HD 18 9733. lPontet dJl2007t HD 209458 ISilva- Valid [2008: 
Tres-1, iRabus et al I [20091 and CoRoT-2, lLanza et al l 12009. 
Especially, iPont et au s study, based on HST data, for the first 
time showed that low-noise transit photometry of exoplanets can 
yield detailed information about surface features of their host 
star. We present a systematic analysis of the spot locations and 
extensions, including their uncertainties, that can be deduced 
from transit lightcurves. 

CoRoT-2a (GSC 00465-01282) is a n apparently young 
late G-dwarf star (lAlonso et al.1 120081 AL08 in the fol- 
lowing); it is unsually active and intrinsically variable 



among the presently known pla netary host stars. Its rota- 
tion period of P, = 4.52 + 0.14 d (lLanza et al.ll2009l) amounts 
to less than three times the planetary orbit period P or b\ 
13 • P or b = 5.004 • P„. Densely sampled high-precision photom- 
etry in combination with spectroscopic measurements make 
the trans it geometry of Co RoT-2b exceptionally well known 
(AL08, Bouch v et al.ll2008l) : The planet is on a nearly circu- 
lar orbit (e = 0.003 + 0.003) which is seen close to edge-on 
(z = 87.84 + 0.10°) and approximately oriented perpendicular to 
the stellar rotation axis (deviating by A = 7.2 + 4.5°). 

These attributes turn the planetary disk into an extremely 
well-defined probe which periodically scans a band on the stellar 
surface covering 20° in latitude. The orbital period of CoRoT-2b 
translates into an orbital angular velocity of 0.002 deg/s. As a 
result, close to the center of the stellar disk, the planet moves 
* 0.1 deg in longitude during a 32 s time bin of the CoRoT data. 
By order of magnitude, this value then determines the attainable 
surface resolution of the CoRoT observations of CoRoT-2. 

2. Observations and data analysis 
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The CoRoT dAuvergne et al.L 120091) lightcurve of CoRoT-2 ob- 
served from 2007 May 16 to 2007 October 15 continuously sam- 
ples 31 stellar rotations. After five days of observations the satel- 
lite's sampling cadence of CoRoT-2 was switched from 512 s to 
32 s, resulting in a lightcurve that covers 79 planetary transits at 
this high time resolution. 

Our analysis started out from the lightcurve data as delivered 
by the CoRoT pipeline. The pipeline sorts out defective data (e.g. 
due to crossings of the South Atlantic Anomaly, SAA) and per- 
forms a background subtraction. As a first step we added up all 
three CoRoT photometry channels (red-green-blue) into a single 
lightcurve, because the individual channels are more affected by 
instrumental instabilities than their summed-up signal. 
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Fig. 1. Section of the CoRoT lightcurve of CoRoT-2, purged of outlying values (see Sect. [2J. The transits appear as recurrent dips 
overlying the modulation by the stellar rotation. Transits are numbered from the beginning of the CoRoT observations, our analysis 
focusses on transit no. 56. 



Similar to AL08, we then removed outlier points which show 
a pronounced non-normal distribution. To this end we computed 
the standard deviation <x in several narrow intervals, yielding 
cr x 950. We discarded all points deviating more than 3<x from 
a boxcar-smoothed copy of the lightcurve, thus rejecting nearly 
2% of all points. Finally, also following AL08, we corrected for 
the 5.6% contamination by a neighbouring object; a subinterval 
of the resulting lightcurve is shown in Fig. Q] 

The transits occur during different stellar rotation phases, i.e. 
at different levels and slopes of the lightcurve. In order to com- 
pare the transit lightcurves, it is convenient to normalize them 
to a common level. To this end, we carried out a linear interpo- 
lation of points adjacent to each transit and divided all transit 
lightcurves by their interpolating function. One lightcurve nor- 
malized in this way is shown in Fig. [2] 

We note that this procedure introduces a systematic error that 
depends on the spot coverage of the stellar disk visible during the 
transit. If, e.g., the disk is covered by dark spots not occulted by 
the planet, the transit depth will not yield the true ratio of radii 
R planet I R*- Instead, the stellar radius R* will be underestimated 
relative to the planet, because part of the stellar disk is dark and 
essentially invisible. Given the typical amplitude of CoRoT-2's 
lightcurve of about 4% peak-to-peak, this introduces a compa- 
rable uncertainty for the planetary radius deduced from a sin- 
gle transit. Consequently, this uncertainty also affects our spot 
size estimate of the next section. We will discuss the influence 
of activity-induced lightcurve variations on the determination of 
planetary parameters in a forthcoming paper. 



3. Transit lightcurve modelling 

To determine a spot configuration compatible with the deforma- 
tions of a transit lightcurve, we selected one particular transit 
occuring close to JD 2454335.0 and referred to as "transit 56". 
Rendered in Fig. [2] it shows the most pronounced and isolated 
"bump" of the whole time series of CoRoT-2, suggesting a rela- 
tively narrow spot occulted close to the disk center. Additionally, 
the symmetry of the bump indicates a spot geometry that is 
largely symmetric in longitude. 

3.1. The model 

To model a transit lightcurve we decompose the stellar surface 
into roughly square elements. The integrated stellar flux is com- 
puted by summing up the fluxes of all visible elem ents, respect- 
ing th eir projected area and the limb darkening (e.g. Wolte retall 
2005). To model the planetary transit, all surface elements oc- 
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Fig. 2. Normalized lightcurve during "transit 56", as a func- 
tion of time from the transit center. Gray symbols indicate the 
unbinned measurements, the red line shows them averaged in 
224 sec time bins with lcr-errors indicated. The vertical lines 
delimit the time interval used for our spot model; the lightcurve 
fit resulting from our model is drawn black (model "BN" in 
Table [TJ. The blue line shows the transit model for an unspotted 
star for comparison. The shallow deformation left of the transit 
center is caused by another spot not included in our model. 



culted by the planetary disk at a given phase are removed from 
the sum, thus treating the planet as a circular disk without any 
intrinsic emission. The surface resolution needs to be sufficient 
to analyze the densely sampled movement of the planetary disk 
over the stellar surface. We choose a surface grid with 750 el- 
ements at the equator, yielding 178 868 elements in total and a 
surface resolution of roughly 0.5° near the equator. Also, this 
high resolution is required to compute lightcurve models suffi- 
ciently free of jitter, not exceeding a few 10~ 4 in our models. 

Next, we adjusted the limb darkening to optimize the fit to 
the transit ingress and egress as well as to a tentative lower enve- 
lope of transit 56 and the surrounding transits. We adopt a linear 
limb-darkening law with e = 0.6; we note, however, that this has 
little influence on our determined spot parameters since the spot 
occultation in transit 56 occurs close to the disk center. 

We compute the position of the planetary disk on the stellar 
disk and the corresponding stellar rotation phase using the or- 
bital parameters and planetary size given by AL08, as well as 
a stellar rotation period of P* = 4.5284 d, determine d by a peri- 
odogram analysis and compatible with the value of lLanza et al.l 
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Fig. 3. Fit quality of our transit lightcurve models as a function 
of spot longitude, shown for the models BC (blue stars), BN (red 
plus symbols) and DS (black curve), of Table Q] Our adopted 
limiting value of x 2 = 2.0 is marked by the dotted line. 



(2009), adopting the lightcurve maximum at JD 2454242.67 as 
phase zero point. To this end, we introduce Cartesian coordinates 
in units of the stellar radius R„ with their origin at the center of 
the star and the observer located on the x-axis. The stellar rota- 
tion axis is assumed to coincide with the z-axis, while the axis 
of the planetary orbit is tilted by 90° - i - 2.2° in the xz-plane. 
Both, the stellar rotation and the planetary orbit are adopted as 
right-handed around the z-axis. 

Using these coordinates, z p i and y p i describe the position of 
the center of the planetary disk projected onto the stellar disk 
visible at the stellar rotation phase <p. The following relations 
apply: 



z p i — cos i ■ a/R„ 



(1) 



where i and a describe the inclination and the semimajor axis of 
the planetary orbit, respectively. 
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where <f> cen is the stellar rotation phase at transit center, while 
Ttrans and P* give the transit duration (first to last contact) and the 
stellar rotation period, respectively. Ttrans = 8200 s can be com- 
puted using Eq. (4) of ICharbonneau et all d2006l) . yrv describes 
the lateral planet position at last contact 
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with R p i giving the planetary radius. Note that Eq. [2] is only 
strictly valid for R t <k a; however, the error is of the order 
of (R t /a) 2 and can be neglected for our analysis. In addition, 
Eq. Q] does not include the projected angle between the stellar 
rotation axis and th e planetary orbital axis A. However, given 
iBouchy et alJ (l2008l) 's value of A « 7° and the proximity of the 
modelled spot to the disk center, this introduces an uncertainty 
of less than 2° for the spot latitude. 



3.2. Transit mapping of a starspot 

Our initial tests showed that the given transit lightcurves do 
not significantly constrain the spot shape and contrast. Hence 
we limit our models to circular spots or circle segments. 
Furthermore, we tentatively adopt two values for the spot-to- 
photosphere contrast: As a "dark spot" scenario we choose a spot 
flux o f 30% of the photospheric flux. Based on Bouch v"et alJ 
(2008)'s effective temperature of CoRoT-2a of 5625 K and 
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Fig. 4. Fit quality^ 2 of our transit lightcurve models as a func- 
tion of the spot radius and its center colatitude. The upper panel 
describes the "bright spot" scenario (models BC, BN and BS of 
Table [TJ, while the lower panel applies to the "dark spot" so- 
lutions (models DN and DS). Subsequent contours mark levels 
of^f 2 = 1,2,5, 10 and 20; the adopted limiting contour,^- 2 = 2, 
is drawn bold red. Because of the near symmetry of the x 2 ' 
contours, only colatitudes in the upper half of the transit band 
are shown. See text for discussion. 



Planck's law at a wavelength of 6000 A, this corresponds to a 
spot approximately 1200 K cooler than the photosphere. This is 
motivated by sp ot temperatures found for other highly active G 
and K stars (e.g. lStrassmeier & Ricdfl998llO'Nea1 1 et alj|2004l) . 
As a "bright spot" we adopt a value of 75% of the photospheric 
flux. This roughly represents the average c ontrast of large spot 
groups on the Sun at visible wavelengths (Wa lton et alJ 12003, 
lAlbregtsen et alJl984l) . Our two spot contra st scenarios are co m- 
parable to those adopted for CoRoT-2a bv lLanza et al.l (|2009). 

Defining the spot contrast and shape in this way, three free 
parameters describe a given spot, namely the spot radius r as 
well as its central longitude <p and colatitude 6. However, as il- 
lustrated by Fig. [3] the longitude is closely confined by the transit 
lightcurve, irrespective of the spot contrast and colatitude. This 
is due to the well defined maximum of the transit bump. Thus 
only two undetermined parameters remain to optimize the fit of 
transit 56's lightcurve: r and 0. We calculated model lightcurves 
for grids in these two parameters, the resulting goodness-of-fit 
X 2 is shown in Fig. [4] 

To calculate^ 2 we rebinned the lightcurve into 224 sec bins, 
and estimated the resulting errors crj assuming Gaussian er- 



ror propagation: x 2 - l/(N - M) ■ Zl =l (Fobs 



-F 



model. 



Here, M = 3 is the number of free model parameters, namely the 
spot longitude, colatitude and radius. N = 8 gives the number of 
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Table 1. Parameters of characteristic spot solutions dis- 
cussed in the text, longitudes and colatitudes are given for 
the spot center. 



Model" 


Long. 


Colat. 


Radius 


Flux" 


x 1 


Area c 


BC 


216.4° 


75.0° 


7.8° 


0.75 


0.8 


0.45% 


BN 


216.5° 


70.0° 


9.5° 


0.75 


1.4 


0.55% 


BS 


216.2° 


80.0° 


8.5° 


0.75 


1.2 


0.47% 


DN 


216.7° 


71.0° 


4.8° 


0.3 


1.7 


0.18% 


DS 


216.2° 


81.0° 


4.8° 


0.3 


1.9 


0.18% 


DEQ 


216.3° 


94.0° 


15.3° 


0.3 


1.6 


0.72% 


" BC, BN and BS stand for 'bright central', 'bri 


ght north' and 



'bright south', respectively. They describe spots with colatitudes 
close to the center of the planetary disk. DN and DS stand for 
'north' and 'south' dark spot solutions, respectively; DEQ repre- 
sents a "dark" spot centered below the equator. 

* Relative to the photosphere 

c Fraction of total stellar surface 



time bins used to constrain the modelled spot, see the interval 
limits indicated in Fig. [2] 

Reducing^ 2 in this way is non-unique because of the signif- 
icant correlation of the spot radius and colatitude, illustrated by 
the slanted contours of Fig. H Tentatively we choose a limit of 
X 1 < 2.0 for models judged as compatible with the observations. 
As an example, Fig.[2]shows a solution yielding^ 2 = 1.4. 

3.3. Results 

As the ^-contours in the upper panel of Fig. H] show, for the 
"bright spot" scenario, the central spot colatitude is confined 
to 6 = 75 ± 6°, i.e. inside the stellar surface belt transited by the 
planet. The spot radii are slightly smaller than, or comparable to, 
the size of the planetary disk: r = 1.1° ... 10.6°. This scenario is 
illustrated by Fig.[3]and the exemplary solutions BN, BC and BS 
("bright north, central and south ") of TableQ] 

For the "dark spot" scenario, on the other hand, only spot 
centers away from the center of the transit belt yield proper fits to 
the transit lightcurve: 6 < 12° and 8 > 18°. This is illustrated by 
the lower panel of Fig. |4]which also shows that the resulting spot 
radii are smaller than in the "bright spot" case (r = 4.6° . . . 7.8°). 

Also spots with centers outside the transit band yield feasi- 
ble solutions. An example is illustrated by the red arc in Fig. [5] 
showing solution DEQ of Table Q] As illustrated by the figure, 
concerning area and longitude extent they do not differ signif- 
icantly from solutions with centers inside the band. We do not 
discuss them further since their radii do not describe the tran- 
sited extension and area of the spot appropriately. 

4. Discussion 

The shape of the transit lightcurves of CoRoT-2 exhibits highly 
significant variations between different transits. This indicates a 
ubiquitous presence of starspots in the surface regions transited 
by th e planet. Furthermo re, as illustrated by Fig.Q]and the analy- 
sis of lLanza et alj (|2009) the overall lightcurve of CoRoT-2 con- 
tinuously changes in amplitude during the complete CoRoT ob- 
servations. This shows that the stellar surface of CoRoT-2a per- 
sistently evolves on timescales shorter than its rotation period. 
Such a fast spot evolution is interesting physically and makes 
CoRoT-2a a favourable case for the study of stellar activity, po- 
tentially a landmark system. 
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Fig. 5. A starspot on CoRoT-2a occulted by the planet during 
transit 56, as reconstructed by transit mapping. The black and 
gray circles represent the planetary disk and the spot for our 
"bright spot" scenario, i.e. adopting a spot flux of 75% relative 
to the photosphere. The purple arcs illustrate the northern- and 
southernmost solutions for this spot flux. The red arc illustrates 
a "dark spot" solution (DEQ in Table [T]), see text for discussion. 

We concentrate our analysis on a single planetary transit 
whose lightcurve shows a pronounced and isolated bump close 
to the transit center. Assuming this bump is caused by a circu- 
lar starspot, we determine parameter ranges for this spot which 
reproduce the observed transit lig htcurve. 

While, similar to iPont et all s analysis of HD 189733, the 
spot contrast is only weakly constrained, the spot longitude and 
radius are closely confined by the transit lightcurve. The spot 
thus reconstructed on CoRoT-2b is comparable in extent with 
large spot grou ps on the Sun which cover up to about 1% of the 
solar surface (Baumann & SolankiLl2^0llNormanLl2005l) . 

Given the nearly 90° inclination of CoRoT-2b's orbit, the 
transit-covered belt on its surface lies close to the equator. Our 
analysis proves that CoRoT-2a, like the Sun, exhibits spots in 
this region. Such well-constrained latitude measurements near 
the equator are difficult or impossible with other surface recon- 
struction methods like Doppler imaging. Using long-term tran- 
sit observations, this may allow to study activity cycles analog 
to the solar butterfly diagram. Also, concer ning possible indica- 
tions of differential rotation on CoRoT-2a dLanza et all [2009), 
the spots in the transit-covered belt could supply additional in- 
formation. 

Solar umbrae have diameters up to about 10 Mm, corre- 
sponding to one degree in heliographic coordinates; penumbrae 
reach approximately twice this size. Our study indicates that a 
surface resolution of potentially better than one degree can be 
achieved for host stars of eclipsing planets when applying transit 
mapping to low-noise and fast-sampled lightcurves. Thus, such 
observations offer a new opportunity to study "solar-like" sur- 
face structures on other stars, they may for example allow to 
measure the umbra/penumbra contrasts of their spots. 

Acknowledgements. U.W. acknowledges financial support from DLR, project 
50 OR 0105. 
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